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Abstract 

We analyze fi + pT pair production generated by high-energy electrons emerging from a laser- 
wakefield accelerator. The u + u~ pairs are created in a solid thick high-Z target, following the 
electron accelerating plasma region. Numerical estimates are presented for electron beams obtained 
presently in the LBL TW laser experiment jl| and possible future developments. Reactions induced 
by the secondary bremsstrahlung photons dominate the dimuon production. According to our 
estimates, a 20 pC electron bunch with energy of 1 (10) GeV may create about 200 (6000) muon 
pairs. The produced fi^ can be used in studying various aspects of muon-related physics in table 
top installations. This may be considered as an important step towards the investigation of more 
complicated elementary processes induced by laser driven electrons. 

PACS numbers: 12.20.Ds,13.60.Le,41.75.Jv 
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I. INTRODUCTION 



100 GV/m with 
2 1 which allows, 



The possibility to produce strong electric fields of the order of 10 — 
present laser facilities is a great advantage for laser-wakefield accelerators 
in principle, to construct compact accelerating devices for particle and nuclear physics. The 
successful production of high-quality electron beams in such laser-driven accelerators with 
electron energies of the order of 1 GeV has been reported recently [l, 3, 4|. Electron beams 
with energies exceeding 1 GeV are interesting for many applications in particle physics, such 
as investigating the properties and production mechanisms of vector and exotic scalar mesons 
in photo/electroproduction 

A3 

, excitation of baryon resonances with the aim of studying 
their properties and search for missing resonances pj, strangeness photoproduction 8] etc. 
Another subject is related to the neutrino physics. For example, neutrino oscillations need 
two types of neutrinos, at least. They may be obtained in muon decays /i + — > e + + u e + 
and pT — > e~ + v e + v IJL , where muons and electrons neutrinos (or antineutrinos) are produced 
in equal parts. Therefore, it is interesting to estimate whether the high-energy laser-driven 
electrons can produce a sizeable amount of muon pairs for future applications. 

Recall that the present generation of the high-energy neutrino beams are made by decays 
of charged pions and kaons in flight in a long decay channel. In this case, the neutrino 
(antineutrino) beams mostly consists of with rather small admixture of electron 

neutrinos. The idea to use a storage ring of muons as a source of high-energy and high- 
intensity neutrino beam has been discussed in Refs. (si, Q . It is assumed that the muons 
are produced in two-body decays 7r + — > p + + and K + — > p + + and then stored in a ring 
with subsequent decay into an electron and two neutrinos. The pions (kaons) in turn, are 
produced in proton-proton and (or) in proton-nucleus collisions with high-intensity proton 



beams 



ll] . Together with neutrino oscillation, such high- intensity muon sources may be 



used in studying other fundamental problems of lepton physics, say the search for lepton 
flavor violation [12j and the measurement of the muon's anomalous magnetic moment 13]. 
The idea to use a laser-driven proton beam for investigating different aspects of the neutrino 



oscillations was discussed for the first time in Ref. 

The aim of present paper is to analyze the possibility of muon pair creation in the 
interaction of high-energy laser driven electrons within a heavy (high-Z) target in a table top 
configuration. The electromagnetic sources of the yU + /i~ (dimuon) production are described 
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by the following elementary processes 

1 + A 
e + A 



e' + A + fi+fi" . 



(1) 
(2) 



In the first case (1) the dimuons are produced in the interaction of real (bremsstrahlung) 
photons within the electric field of the high-Z target nuclei. This is an analog of well 
known Bethe-Heitler process of the electron-positron production. In the second case (2), 
the dimuons are produced in the interaction of high-energy electrons with nuclei (so called 
trident process). These two reactions are depicted in Fig. 1. In some sense these reactions 



1-1+ 



n- 2 



(a) 



(b) 



FIG. 1: Diagrammatic representation of dimuon production processes in electromagnetic in- 
teractions, (a) Bethe-Heitler process 7 + A — > A + (b) Electron induced reaction 
e + A -> e + A + fi + ^~. 

are similar to the dielectron production by laser-driven relativistic electrons at few dozen 
MeV analyzed by Nakashima and Takabe 15|] and other authors (see, for example, Ref. 16]). 
Therefore, our investigation may be considered as a continuation of these studies, but with 
the focus on dimuon production. For other avenues towards laser-generated muon pair 
production cf . Refs. [o, Q • Dimuon production in electron-positron-photon plasmas was 



estimated in 



19] 



Our paper is organized as follows. In Sect. II, we recall the relevant elementary elec- 
tromagnetic processes. Estimates of the dimuon yield are presented in Sect. Ill, where 
the energy loss and attenuation of the primary electron beam and the secondary photon 
bremsstrahlung spectrum are accounted for. Hadronic processes are briefly discussed in 
Sect. IV. Our conclusions can be found in Sect. V. 



3 



II. ELEMENTARY CROSS SECTIONS 



The differential cross section of the dimuon production in the reaction Eq. ([I]) reads 



da^ A = M | T ^| 2 yj M* +I1 _ - Ay? dM^+y- d cos e q d% + , (3) 

where k and q are the spatial parts of the incoming photon momentum and the total mo- 
mentum of the outgoing muon pair q = + p„- with the invariant mass of the muon pair 
M 2 +At _ = q 2 ; fi denotes the muon mass, 6 q is the polar angle of the direction of flight of 
dimuon, is the solid angle of the direction of the velocity of the fi + meson in the rest 
frame of the dimuon; Z and a stand for the nuclear charge and the fine structure constant 
(a = 1/137), respectively; T lA is the invariant amplitude and averaging over photon polar- 
ization and summation over spin states is provided in |T 7j4 | 2 . The quantization axis z is 
chosen along velocity of incoming photon, and the production plane is defined by the vectors 
k and q. The differential cross section depends on the initial energy and four kinematical 
variables, which define the final state. 

In case of the e + A — > e' + A + fi + fi~ reaction Eq. ([2]) we have relevant two sub-processes. 
One corresponds to the electron scattering off the nucleus. In the second case, one of the 
outgoing muons interacts with the atomic nucleus. These two cases are marked in Fig. 1 (b) 
as " 1" and " 2" , respectively. The total amplitude is the coherent sum of all four amplitudes. 
The differential cross section for this reaction depends on seven kinematical variables, which 
define the final state of the outgoing leptons and recoil nucleus, and it has the form 

'2^,4 



eA _ Z a |p e /| , A|2 



= 7TT \T eA \ 2 ^M^ +ii _ - 4/i 2 \q\dE q dM^- dcos6 e , dQ q dty, + , (4) 



where p e and p e ' are three momenta of the incoming and outgoing electrons, respectively, 8, 



is the polar angle of the direction of flight of the outgoing electron, E q — q 2 + M 2 +fi _ is the 



energy of the muon pair; T eA is the invariant amplitude and averaging over initial electron 
spin projections and summation over the final fermion spin states is understood in |T eA | 2 ; the 
quantization axis (z) is along velocity of the incoming electron, and y = [p e x p e /]/|p e ||p e /|. 

The invariant amplitudes for reactions in Eqs. ^ and ([2]) are calculated in the lowest 
order of perturbation theory analog to the electron-positron production in 7 A and eA 
interactions [20]. The amplitude for the 7 A reaction reads 
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where u and v are the Dirac spinors of the outgoing fi~ and /i + muons, respectively, e 7 is the 
polarization vector of incoming photons, g\ = k — p^+ , g 2 = k — p^- , and k 7 = k — q. The 
notation f means the four product 7ojo Q , where 7 Q (a = 0, 1, 2, 3) are Dirac's 7 matrices. 
Correspondingly, for the e A reaction we have 
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(6) 
(7) 



where q' = p e - p e ,, Pl = p e - q, p 2 = p e , + q, p 3 = q' - p +, p 4 = </ - p _, and K e = q' - q. 



The effects of nuclear size and higher-order Coulomb corrections are small 2l|; the effects 



of the atomic electron screening 22j, |23j are neglected here as the production processes 



happen in the high-field region near the target nuclei. The corresponding cross sections 
Eqs. ([3]) and (j3j) are evaluated numerically without any additional approximation. 
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FIG. 2: Differential cross sections of dimuon production in eA (left panel) and 7 A (right panel) 



interactions as a function of the invariant dimuon mass M, 
correspond to initial energies of 1 and 10 GeV, respectively. 



The solid and dashed curves 



The invariant-mass distributions for 7 A and eA reactions are shown in Fig. [3 In this 
case, the corresponding cross sections are calculated as a function of the invariant mass 
M^+p- integrating over the other variables in Eqs. (J3J) and (HI). All calculations are done for 
a gold target. The initial electron energies are chosen to be 1 and 10 GeV. For 1 GeV we 
follow the laser accelerator conditions of Ref. The case of 10 GeV may be considered as 



prediction 



now 



24 



25 



or 



'uture high-energy table top electron accelerators being under consideration 
26] . For simplicity, in case of the 7 A reaction, the elementary cross sections are 
calculated at E 1 = 1 and 10 GeV. Below, for the estimate of the dimuon yield, we fold the 
elementary cross section with the bremsstrahlung photon distribution for all kinematically 
allowed photon energies. 
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FIG. 3: Differential cross sections of dimuon production in e A (left panel) and 7 A (right panel) 
reactions as a function of polar angle 6 q . The solid and dashed curves correspond to initial energies 
1 and 10 GeV, respectively. 

Fig. [2] demonstrates that the cross section of the electron induced reaction is smaller by 
two order of magnitude than the corresponding cross section in the photon induced reaction. 
This is mainly because of an additional vertex with additional factor a. Both distributions 
exhibit a maximum slightly above ~ 2/x. This means that the relative kinetic energy of the 
two muons in a pair is small, i.e. they flight close to each other. 
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FIG. 4: The differential cross section for electron induced muon pair production reaction as a 
function of dimuon energy E q . The solid and dashed curves correspond to initial energies 1 and 10 
GeV, respectively. 
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The differential cross sections as a function of the polar angle of the direction of flight of 
the muon pair 9 q integrated over the other variables are shown in Fig. [3j One can see that 
the cross sections are peaked at the forward direction, especially for E en = 10 GeV. The 
knowledge of the spatial structure of the outgoing muon flux is important for the design of 
the devices for muon sources. 

The differential cross section of the reaction e + A — > e' + + A as a function of the 

dimuon energy E q is presented in Fig. HI One can see a wide spread of this distribution with 
a maximum around 1 GeV for initial electron energy of 10 GeV, while for 1 GeV electrons 
the dimuon energy displays a plateau above 0.6 GeV till the kinematical limit. Since the 
kinetic energy of the relative motion of /i + and /i~ within a pair is small {E^± ~ E q /2), the 
energy distribution of the individual muons may be approximated as da/dE^+ ~ da jdE^- ~ 
da/dE q \ Eii± =E q /2- 
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FIG. 5: The total cross section of the dimuon production in 7 A (solid curve) and eA (dashed 
curve) reactions as a function of initial energy. The dotted curve corresponds to the high-energy 
approximation Eq. ([8]). 



The total cross sections of the dimuon production as a function of the initial energy 
are shown in Fig. One can see a steep increase of the cross sections with energy. We 
note that in the ultra-relativistic case with E^ ^> fi the total cross section of the reaction 
^A — > ji + \T A may be described by 



28 2 m2 / 2£ 7 109\ 



(8) 



where r§ = a/fi is the classical muon radius. The corresponding cross section is shown in 
Fig.[5]by the dotted curve showing that the ultra-relativistic behavior sets in at E 1 > 3 GeV. 



The total cross sections for 7 A and eA reactions reach 300 fib and 3 fib, respectively, 
at i? 7 (e) — 2 GeV. Similarly to the differential cross sections, the difference between the two 
reactions is about two orders of magnitude, mainly due to an additional power of a in the 
e A cross section. 



III. DIMUON YIELDS 

Using the above elementary cross sections one can estimate the dimuon yield for given 
electron beam and target properties. For the former ones we use the conditions of electron 
beams as reported for the laser- wakefield accelerator in Ref . l| . The electron energy is about 
0.5 — 1 GeV and the electron flux is about 20 pC which corresponds to N$ ~ 1.248 x 10 8 
electrons in a bunch. In our estimates we assume the same flux for electron energies up to 
10 GeV. We consider a gold target with thickness of L = 0.1 — 1 cm. 

Strictly speaking, the particle production in interactions of the high-energy electrons 
with heavy target nuclei must be evaluated by transport-kinetic models (see, for example, 
, ll6|). However, for a first qualitative estimate of the dimuon yield one can use an 



Refs. 
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analytical approach, similar to that developed in Ref. [2j 

Consider first the dimuon yield which stems from elementary eA — > e' + fi + fi~ + A 
reactions. It may be expressed as 

L 

dN ^+^- = ^A1A J dlN e (l)da eA ^ e,fl+f *~ A (E e (l)), (9) 


where A is the atomic weight, Na is Avogadro's number, pa denotes the target density, 
d(J eA-*e' ti+ n a j g ^ e elementary cross section of the dimuon production discussed above. For 
the sake of simplicity, we neglect the energy spread in the bunch taking the electron energy 
at the central positions of the energy distribution. This seems to be reasonable, because the 
energy spreading reported in l] is less than 100 MeV at E e ~ 1 GeV. 

Propagating through the target material the electron beam loses its intensity and energy. 
The energy loss is described by 

dE e _ NaPaoZ^I ( 183 2\ 

with Tq = a/m e ~ 2.18 x 10~ 13 cm. For the gold target with A = 197 and Za = 79 this 
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leads to the formula 

E e (l)~E° e exv(-l/l ) (11) 

with Iq ~ 0.513 cm. The electron beam absorption in the target may be accounted for by a 
linear dependence 

iV e (/)^iV e (l-//L max (E e )) , (12) 

where Nq is the initial beam intensity and L max stands for the maximum distance traveled by 
the high-energy electrons. For the gold target one can use L max ~ 1.44 cm at E e = 0.5 GeV 



which may be extrapolated as L max (E) ~ 0.23 E ' 3 cm with E in GeV 28|. Both, the 
decrease of the beam intensity and the energy loss lead to a decrease of the effectiveness 
of the dimuon production in the target material. However, the effect of the energy loss is 
significantly greater than the effect of intensity depletion since the scale parameter l in 
Eq. (II ip is much smaller than L max . 

The yield of dimuons, produced by secondary real photons created at the distance I 
behind the target front side, reads 

L 

dN ^- = ^AfA f ^ I dE ^^M± {L _ /} da jA^-A( E j ? (13) 

A J J dE^ 



o 

where da' yA ~^ fl+fJ ' A is the above cross section of the dimuon production in the elementary 
7 A interaction, and dN^(l) / l dE 1 denotes the distribution of bremsstrahlung photons with 
energy E 1 at distance / 

dN^l) _ N APA ATein da^E e (l),E y ) 
dE 1 ~ A [) dE y ' 1 > 

where N e (l) is defined by Eq. (1121) and the current electron energy E e (l) from Eq. (II ip is to 

be used. For the angular-integrated bremsstrahlung cross section daJE e , E 1 )/dE 1 we adopt 

a parametrization motivated by the analytical expression of Ref. [20j 

fjg = ia < F ^ {(El + El - \E^) in + E -f-) , (15) 

where E e , = E e - E 1 and F{E e ,E 1 ) = 0.91 tanh(£ e ,/0.02£ e ). In Fig. E we show the 
bremsstrahlung cross section for two energies together with the compilation of Ref. [29 ]. 

Results of calculations of the dimuon yields iV M+M predicted by Eqs. (Q and (|T3|) as a 
function of the primary electron energy E® are exhibited in Fig. [71 The total number of 
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FIG. 6: The scaled bremsstrahlung cross section (E~ f /Z 2 )da~ f (E e ,E~ f ) as a function of E^/E e for 
two electron energies 1 GeV (solid curve) and 10 GeV (dashed curve). The symbols represent the 



data compilation of Ref. 



29]. 



dimuons produced in e A interactions is about 1 and 60 for E® = 1 and 10 GeV, respectively 
for thick target. The yield increases with the thickness of the target. Some saturation sets 
in at L > L sat ~ 0.5 cm. This means that the dimuons are essentially produced in a narrow 
region of the target mostly because of the large energy loss. It seems to be natural, because 
the saturation length L sat . is close to the scale parameter / m Eq. (jPTI) . 

The number of dimuons in 7 A interactions is about 200 and 6000 for E® = 1 and 10 GeV, 
respectively, for target thickness of 1 cm. The dependence of the total yield on the target 
thickness is rather strong due to the large main free path of the bremsstrahlung photons. 
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FIG. 7: The yields of dimuons in e A (left panel) and 7 A (right panel) interactions as a function 
of primary electron energy E® , calculated by Eqs. © and (fT3|) . respectively. 



The dimuon yield in 7 A interactions considerably exceeds the corresponding yield in e A 
interactions. This excess strongly depends on the target thickness. Therefore, in Fig. [8] 
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we present the total dimuon yield in interactions of relativistic electrons with a gold target 
which is a sum of the two above contributions as a function of the primary electron energy 
Eg and target thickness L. For thick targets it practically coincides with the result shown 
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FIG. 8: The total yield of dimuons in interactions of relativistic electrons with the gold target as 
a function of the primary electron energy E® and target thickness L. 

in Fig. [7] (right panel). 




IV. BRIEF COMMENT ON HADRONIC PROCESSES 

Finally, let us consider the hadronic sources of muons in interactions of the laser-wakefield 
accelerated electrons with the thick target. The main source of muons here is the photo 
and/or electro production and subsequent decay of 7r =l= mesons in the elementary processes 

7 + p — > 7r + + n, 7 + n^7r~+p, (16) 
7 + iV->7r + + 7r~ + iV, (17) 

where 7 stands for real bremsstrahlung photons (photoproduction) or virtual photons (elec- 
troproduction). The production of K mesons by 1 GeV electrons is strongly suppressed 
because of energy arguments, however, at higher energies they contribute to the muon pro- 
duction similarly to the pions. Charged pions decay via 

7T + —> fi + + Up, 7I~ -> p~ + Up (18) 

with a probability of 99.99%. The cross sections of the reactions (TIE]) and (TIT1) for E 1 ~ 
1 GeV are around 100 pb (see, for example, [so|). Pion photoproduction in 7 A reactions is 
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proportional to A 2 / 3 . This means that the cross section of the muon production at heavy 
target nuclei by 7 A — ► n^X — > yU ± +z/ M (z/ At )X reactions exceeds the dimuon photoproduction, 
considered in previous section, by a factor of 20 — 30. Taking into account the large pion 
absorption cross section at heavy target nuclei by strong interaction processes this means 
that the produced pions will be essentially absorbed in a thick target. The target must be 
thin enough to get a well collimated muon beam. This condition reduces the muon yield 
considerably. For example, for a dozen /zm thick gold target the expected yield of ^ muons 
is less than 1 event for the above considered electron bunch. Therefore, the direct production 
of dimuons in electromagnetic interactions seems to be a more favorable mechanism. 



V. SUMMARY 



In summary we have considered muon pairs production by GeV electrons, created by a 
laser wakefield accelerator, impinging on a thick high-Z target. We estimated the effective- 
ness of a such a source of muon pairs. For a 1 cm thick gold target, 1.25 x 10 8 electrons 
in a 20 pC bunch with energy of 1 (10) GeV in the initial state produce about 2 x 10 2 
(6 x 10 3 ) dimuons with pair energies centered at 1 GeV. To get 10 6 dimuons from the muon 
factory one needs 10 10 — 10 11 primary electrons in a bunch. Such intensities with power of 
100 J seem to be quite realistic in near future, requiring ultra-high intensity laser pulses with 
efficient acceleration mechanisms. Cooling of the heat load in the target material analog to 



Ref. 



3 if ] may be an option for higher repetition rates. The produced muons, unlike elec- 



trons and hadrons, penetrate the target material without suffering noticeable scattering and 
absorption. Thus the configuration of a laser driven electron accelerator and thick high-Z 
target may serve as an all-optics table top device for muon pair production. The produced 
(jl mesons may be used in studying various aspects of muon and neutrino physics and to be 
considered as an important step towards investigations of more complicated electron induced 
elementary processes. 
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